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Abstract 

The aim of an invisibility device is to guide light around any object 
put inside, being able to hide objects from sight. In this work, we propose 
a novel design of dielectric invisibility media based on negative refraction 
and optical conformal mapping that seems to create perfect invisibility. 
This design has some advantages and more relaxed constraints compared 
with already proposed schemes. In particular, it represents an example 
where the time delay in a dielectric invisibility device is zero. Further- 
more, due to impedance matching of negatively refracting materials, the 
reflection should be close to zero. These findings strongly indicate that 
perfect invisibility with optically isotropic materials is possible. Finally, 
the area of the invisible space is also discussed. 
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1 Introduction 



In Optics, Fermat's principle describes the trajectory of the light as the 
path taken between two points by light waves that can be traversed in the 
least time [TJ |2] . Each dielectric medium is characterized by a refractive 
index n. This index profile integrated along the ray trajectory defines the 
path length. The value of the refractive index indicates a property of a 
material that changes the speed of light, computed as the ratio of the speed 
of light in a vacuum to the speed of light through the material. When 
light travels at an angle between two different materials, their different 
refractive indices determine the angle of transmission (refraction) of the 
ray of light. This relationship is known as Snell's law and illustrates 
that light changes the trajectory to minimize the time for travel between 
two points. A well-known example of the optical effect caused by light 
passing through two media with different refractive indices is a mirage in 
the desert. The difference in air density changes the refractive indices and 
makes it possible to see images from the sky above the sand. This is a 
consequence of the Fermat's Principle, because the light rays follow the 
shortest path defined by the lowest refractive profiles. 

Our invisibility problem has strong links with these principles of Op- 
tics. Recently, a few models have been proposed to create a perfect illusion 
of invisibility. The main idea is to design a dielectric medium composed of 
a specific refractive index that literally bends and guides the light around 
a desired object. As the device itself should be invisible, an external ob- 
server would not see the object. The optical effect will be equivalent to 
observe the light rays propagating across empty space, because objects in 
the background will be visible. Until now, a perfect invisibility effect based 
on isotropic media has not been achieved and the ultimate wave nature of 
light has been argued as the main reason. This issue has been extensively 
debated in several works |4] [5] |6] . However, using anisotropic media 
these distortions could be reduced to zero, in principle 17 . Furthermore, 
a recent experimental demonstration of a cloaking device based on arti- 
ficially structured metamaterials [8] has succeed in decreasing scattered 
waves emitted from the hidden object. The object remains invisible over 
a narrow band of microwave frequencies. More recently, a new design 
of a non-magnetic optical cloak that reduces the scattered waves, and 
extends the operative range of frequencies to the optical spectrum, was 
proposed and investigated using computer simulations [9]. Implications 
of these research lines, future trends and remaining problems on the field 
are commented and summarized in |10| . 

The discovery of new materials (also known as metamaterials) is open- 
ing up exciting avenues for new electromagnetic applications in diverse 
areas, such as medical imaging, microwaves, lenses, radars, defense and 
telecommunications. These new features cannot be derived using mate- 
rials found in nature. Metamaterials then represent artificial media with 
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unusual electromagnetic properties. These material properties emerge as 
a consequence of their complex and often periodic structures rather than 
their chemical composition. Therefore, engineered materials with pre- 
designed periodic nanostructures could lead to a large range of refractive 
indices with a rich and great variety of new electromagnetic properties. 
This is also true for the recently discovered metamaterials with negative 
refractive index [llj . A striking property of materials with negative in- 
dex is that the rays will be refracted on the same side of the normal on 
entering the material according to the SnelPs law. The connections be- 
tween invisibility and negative refraction were described in [12] using the 
formalism of general relativity. 

In this work, we extend the original idea of one of us based 
on optimal conformal mapping to create dielectric invisibility devices. In 
this context, we propose a new design for a cloaking device made of a 
metamaterial with negative refractive index. Our scheme reduces the time 
delay caused by the device to zero and relaxes the constraints imposed by 
the requirement for bounded orbits in the Riemann sheet. By following [5] 
|6] , a dielectric medium conformally maps a physical space z onto Riemann 
sheets given by an analytic function w(z). As a feature of conformal maps, 
the angles between coordinate lines are conserved. By using, for example, 
the well-known Joukowski transformation 13 , 14 , the papers exploit 
a mapping between the physical space described by the complex field z 
and the analytic function w(z) that represents the mathematical space 
composed of two Riemann sheets and a branch cut which connects both 
sheets. As a result any object located at the center of physical space 
cannot be detected by an external observer. In other words, it seems 
hidden inside the Riemann sheets. 

The most important distortions of invisibility are caused by reflections 
and time delays [6]. In this work, we investigate how to relax the con- 
straints for designing and manufacturing invisibility devices. In previous 
works, refractive-index profiles on the interior Riemann sheet should guide 
the rays around the branch points and should define close loops or tra- 
jectories. Therefore, according to classical dynamics [2], the number of 
potential or index profiles is drastically reduced to a few ones: harmonic- 
oscillator, Kepler profile and Maxwell's fish eye. In contrast, our proposed 
design is able to generate a large variety of bounded orbits for rays us- 
ing striking properties of negative refractive indices. More importantly, 
the proposed device based on negatively refracting material represents an 
example where the time delay in a dielectric invisibility device is zero. 
Furthermore, due to impedance matching of negatively refracting materi- 
als the reflection should be close to zero. These findings strongly indicate 
that perfect invisibility with optically isotropic materials is possible. Fi- 
nally, the area of the invisible space is also discussed. 



3 



2 Mathematical Formulation 



In this section, we follow refs. [5] \G\ and describe the rays of light based 
on the Hamilton's analogy I I . between the trajectory of rays in media and 
the motion of particles governed by classical mechanics principles. Let us 
assume that the refractive index n(x) does not vary too much comparing 
with the wavelength of light. 

Then, the electromagnetic wave in media can be expressed using the 
Helmholtz equation: 



(A+^)^ = 0. (1) 

c 

Next, we introduce the effective time r measured in spacial units 

dr = ±dt, (2) 
n z 

it is possible to derive the trajectory of rays as follows [2) [5]: 

-Vn 2 . (3) 



(fx 1_ a 



dr 2 2 

This equation resembles the Newton's second law equation for light rays. 
Integrated ^ by x, we obtain 

1 , dx . 2 n 2 . . 

The corresponding energy E and potential U follows the relationship [TJ[5]: 

n 2 

U-E = -\. (5) 

Here, we can show that this constant in the right hand side of @ vanishes 
as follows: 

Therefore, from (O we obtain 

\f T \=n. (7) 

In what follows, we use the Newton's second law equation for light rays 
© for designing invisible devices. 



3 Conformal mapping 

In this section, we review the conformal mapping introduced by [6] 
in the context of invisibility. This mapping is used to design and create 
the invisible space inside the material. The property of this conformal 
mapping and the associated Riemann surfaces plays a crucial role to design 
the cloaking devices. 
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3.1 Complex plane 

We assume that the medium is uniform along the z-direction. Then, we 
introduce the complex coordinate z = x+iy and z = x — iy. The derivative 
of complex coordinate is given by 

d~z-2 { d^ + lW' (8) 

dz 2 y dx idy>' K > 

Next, using complex coordinates, Helmholtz equation |T} is transformed 
into 

Let us now consider the analytical function u> = w(z), which is con- 
formal mapping. Then, by using 

dz dz dw 

d_ = dwd_ 

dz dz dw ' 

we find that, by the conformal mapping w = w(z), the Laplacian operator 
is transformed as follows 

d 2 ,3to, 2 d 2 



(13) 



dzdz dz dzdz 
When the refractive index is transformed as 

n z = n w \ — | (14) 

by conformal mapping, the Helmholtz equation (|10p is invariant. Here, 
n z and n w are the refractive index of z and w complex plane, respectively. 
In what follows, we use this invariance to design the invisible device. 

3.2 Conformal mapping and Riemann surface 

We map the physical space z into Riemann surface w by using the follow- 
ing conformal map: 

w = z+— . (15) 

z 

The inverse function of (|15[1 is given by 

z = i(«>± Vw 2 -4a 2 ). (16) 

Here, w-space is a Riemann surface which consists of two sheets. By the 
conformal map 1)150 . the outside of the circle |z| > a is mapped into the 
first Riemann sheet in w-space, and the inside of the circle zj < a is 
mapped into the second Riemann sheet in ui-space. The branch cut is 
given by —2a < Re(w>) < 2a, lm(w) — in m-space, which corresponds to 
the circle \z\ = a in z-space. 
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3.3 Property of conformal mapping 

When we take \z\ —+ oo, the conformal map (|15p is simplified as follows: 



w ~ z. (17) 

Therefore, in the asymptotical region 2 — > 00, z-space and the first 
Riemann sheet of w-space can be identified with each other. 

We set z as polar cordinate z = re lS and w as cartesian coordinates 
w — u + iv, then the conformal map (|15[1 maps the circle (r = ro) in 
physical z-space into the ellipse: 



(a 2 /r + r ) 2 (a 2 /r -r ) 2 
in tu-space. This property will be used later. 



4 New model of cloaking devices 

4.1 Negative-refraction metamaterials 

Metamaterials are man-made media where electromagnetic waves do not 
behave as usually expected. As a consequence, new electromagnetic prop- 
erties and effects may emerge and be observed. This term is particularly 
used when the material has properties not found in naturally-formed sub- 
stances. Therefore, engineered materials with complex and pre-designed 
nanostructures could lead to a large range of refractive indices with a rich 
and great variety of new electromagnetic properties in many diverse fields 
as microelectronics, medical imaging as well as defense and telecommuni- 
cations [7J[TT]. This is also true for the recently discovered metamaterials 
with negative refractive index. A striking property of materials with nega- 
tive index is that the rays will be refracted on the same side of the normal 
on entering the material according to the Snell's law. Then, it allows to 
use this property to design novel cloaking devices. Here, we propose a 
device with two Riemann sheets. While the first Riemann sheet the ma- 
terial has refractive index n=l, the second Riemann sheet is composed 
of four quadrants with alternative values negative and positive refractive 
index profile n. As a consequence, the existence of quadrants with dif- 
ferent refractive index makes it possible to generate bounded orbits with 
more relaxing constrains than suggested previously [SJ [B] . By using the 
negative refractive index, the trajectory around branch cut points in the 
second riemann sheet can be seen as a mirror- like reflection of light. We 
will explain this in detail in the next section. 

4.2 Refractive index 

We first define the refractive index n w on Riemann suface w defined by 
the conformal mapping (|15[) . We can obtain the refractive index n z on 
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physical space z by using transformation shown in (|14p . In what follows, 
for the matter of convenience, we write just n for both refractive index n z 
and n w , when the meaning is clear from the context. 

The space w is composed by two Riemann sheets. The refractive index 
of the first Riemann sheet in space w is given by 

n w = 1. (19) 

Next, we divided the second Riemann sheet in space w into four regions 
as follows: 

Region 1 (Rl): (x < 2a and y > 0), (20) 
Region 2 (R2): (x > 2a and y > 0), (21) 
Region 3 (R3): (x > 2a and y < 0), (22) 
Region 4 (R4): (x < 2a and y < 0). (23) 

The refractive index on the second Riemann sheet in space w is given by 

n w = ^2E - 2U(x) 



' y/2E - 2A(2a - x) (Rl : x < 2a and y > 0) 
- v /2£ + 2.4(2a-:r) (R2 : x > 2a and y > 0) 

y/2E + 2A(2a - x) (R3 : x > 2a and y < 0) 
-sJIE - 2A(2a - x) (R4 : x < 2a and y < 0) 



(24) 



These equations implicitly define the potential U(x) of the Newton equa- 
tion for the light ray. 

It is worth noticing that the refractive index on the second Riemann 
sheet in space w is completely antisymmetric with respect to the line 
x = 2a and s-axis. 



4.3 Energy constraints 

Since the value inside the square root in the right hand side of (|24p in Rl 
region should be positive (n 2 = 2E — 2U should be positive), we obtain 
the constraint 

2E > 2A(2a - x). (25) 

In order to satisfy the previous condition (|25p for any starting position on 
branch cut —2a < x < 2a, y = 0, we have the constraint 

E > 4Aa. (26) 

Therefore, from (|26[1 . A is ranged between 

< A < — . (27) 
4a 

In what follows, we fix the value E, and consider A as a free parameter 
restricted by (g7J . 



7 



4.4 Trajectory of light rays 

We will compute the trajectory in Rl region. The other three regions can 
be obtained by taking a mirror-like reflection of the Rl region. This is 
the main idea of our approach. 

Next, we solve the Newton equation ((3]l under n w given in (|24p . For 
the Rl area, the trajectory is given by 

x = ^Ar 2 + v x r + x , (28) 

V = VyT. (29) 

Here, v x and v y are the initial velocity at the branch cut in second Rie- 
mann sheet w and the absolute value of velocity is given by 

vl+v 2 y = n =2E -2A{2a~x Q ). (30) 

The novelty is that, using the negative refractive index, we can enclose 
the trajectory in the second Riemann sheet in space w (see Fig. HJ. 
Therefore, we can make the starting point and ending point be the same. 
Moreover, the direction of initial vector and final vector can coincide with 
each other. 

In contrast, [5] |6] proposed a few examples to enclose the trajec- 
tory using bounded orbits according to classical mechanics: Kepler and 
harmonic-oscillator potentials, and Maxwell's fish eye. We think that our 
approach is more flexible and admits a large variety of orbits by using 
the negative index profile features. Thus, it is not limited to the small 
number of bounded orbits described by classical mechanics. 



4.5 The maximum point in the second Riemann 
sheet 

Let n be the parameter when the light arrives to the boundary between 
Rl and R2 area. This can be obtained by solving 

1 2 

-Ar a + v x n + x = 2a. (31) 
The solution of this equation (|31[1 is given by 



-v x + V^l + 2A(2a - x ) 

n = -r • (32) 

From (I28|) . we obtain the minimum appoint of the parabolic curve in the 
721 region as follows 

v 2 

Xmin = ~-^+XQ- (33) 

From (|29[) . the maximum y-point in Rl region is given by 



-Vx + \A! + 2A(2a - x ) 

y m a X = VyT! = Vy . (34) 
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In what follows, we consider the case that the light enters vertically v x = 
for the matter of convenience. Then, (|34|l becomes 



Umax = 2^ - (2a - x ))(2a - x )- (35) 

We will maximize y m ax, by changing the starting point xo- This can 
be done by consider the following two cases (i) and (ii), according to A 
value. 

(i) In case of (g < A < £) 

In case of (g| < A < jM, (|35l) takes the maximum 

»— = ! (for £ <j4< l } (36) 



when the light ray starts at xo = 2a — E/(2a). 
(ii) In case of (0 < A < ^) 

In case of (0 < ^4 < J|), (|35p takes the maximum 



E E 
ymax = 4 J— a - 4a 2 (for < 4 < — ) (37) 

V A OQ 

when the light ray starts at xo = —2a. This starting point is the left hand 
end point of the branch cut. 

In both cases (i) and (ii), with increasing A, y max in space w decreases, 
and the invisible space of z becomes larger. We will discuss this issue in 
more detail in section [7] 



5 Time delay 



Due to the opposite signs in the refraction index described in the second 
Riemann sheet, the proposed media presents by construction an example 
where the time delay in a dielectric invisibility device is zero. In particular, 
it applies to devices based on optically isotropic materials. 

The time delay is conformal invariant and can be computed by 

cT = J n z \dz\ 

i w \dw\. (38) 

Here, we compute the time delay. First, we focus on the Rl region. 
Then, the time delay for Rl region is given by 



cT R1 = / n w \dw\ 



Ci 
n 



= -A 2 t? + Av x t? + {2E-2A(2a-x a ))r 1 , (39) 
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where Ci is the trajectory of the light in Rl region. Inserting (|32|l into 
this equation (|39p . we obtain an explicit form as follows 



2 2 4 „, . , + Jvl + 2A(2a - x ) 2 , , , . 

C T fll = (2£ - - -A(2a - xo)) V ^ + 5". (2a - x). (40) 

We consider the other regions. The time delay on R2 region is the 
same absolute value as on Rl region but with opposite sign, because the 
refractive index on R2 has the opposite sign of that on Rl, 

cTri = / n w \dw\ = — n w \dw\ = — cTr 2 (41) 
J C\ J C2 

where C2 is the trajectory of the light in R2 region. 
Similary, for R3 and R4 region, we obtain, 

cT RZ = f n w \dw\ = - f n w \dw\ = -cT R4 (42) 

where C3 and C3 are the trajectories of the light in R3 and R4 region 
respectively. 

Therefore the total of them is completely cancelled out, and the total 
time delay is exactly zero. 

cT to tai = cTm + cTr2 + cTr3 + cTra 

= (43) 



6 Reflection 

Together with time delay, reflection waves are also a distortion that pre- 
vents to achieve a perfect invisibility device. However, in the proposed 
device, due to impedance matching of negatively refracting materials, the 
reflection should be close zero. This finding, together with the zero time 
delay shown in the above section, strongly indicates that perfect invisibil- 
ity is possible in isotropic media. 

New materials called metamaterials have recently been created. They 
are characterized by a negative refractive index. These materials can 
provide a total refraction phenomena when the wave impedances of the 
two media are matched. As a consequence, there is no reflected wave, and 
the distortions due to reflection are zero. Our proposed invisibility device 
consists of a plane with positive/negative n boundary. In this case, the 
reflection would be close to zero by impendance matching. To be precise, 
given two materials with a defined border, impedance matching means 
that Zi n =Z ou t i.e., output impedance of a source of waves is equal to the 
input impedance of the waves, and there is no reflected wave |15| . 

Furthermore, it is worth noticing that our case does not contradict 
Nachman's theorem about the impossibility of perfect invisibility in stan- 
dard isotropic materials, because these negatively refracting materials ex- 
ploit the polarized nature of electromagnetic waves. It means that these 
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materials are described in terms of the permittivity e and the permeability 
fi. 

Let us illustrate the above described reflectionless property of our pro- 
posed device using an incident plane wave solution of Maxwell's equation 
propagating through a planar dielectric interface. 

Let us first consider the reflection coefficients that can be obtained 
from Maxwell eqs. and the boundary conditions at dielectric interfaces 
(see [15] for details). In case that E is vertical to the incoming plane the 
coefficient reads as 

Er _ y^cos^y^/l-^sin^ 
t? r~r i ' ^ > 



where 4> denotes the angle of incidence. In case that E is parallel to the 
incoming plane the coefficient is given by 



Er _ \/f co ^-y^\A-(^) 2 ^ (45) 



For a general angle of E, one should consider linear combination of 
these two cases. There are three limiting cases (1), (2) and (3): 

(1) e = —e and fj, = —fjf case. In this case, there is negative refraction 
and no reflection at all for a general angle of incidence and polarization: 

j = 0. (46) 

We remark that this holds for any polarization, by considering the linear 
combination of the two cases (i) E is parallel to incoming plane and (ii) 
E is vertical to incoming plane. 

(2) Normal incidence case. In this case, the coefficient reads as 



Er _ V e V_£ _ Z 2 — Zl 
Ei l~E _|_ ft[ %2 + Zl 



(3) Impedance- matched case. When Z — Z' , the coefficient reads as 
rp cos A — , 1 — (-^) 2 sin 2 6 

f = V y " ■ (48) 

* cos0+ Jl - (^) 2 sin 2 <^ 

In our problem, at the boundary of Rl and R2, we can take the fol- 
lowing form, i.e., case(l): 

em = — ej?2 (at the boundary between Rl and R2) (49) 
Uri — —/J.R2 (at the boundary between Rl and R2) (50) 

where em and €R2 are the permittivity of Rl and R2 region respectively, 
/Xfli and hr2 are the permeability of Rl and R2 region respectively. This 
gives the property 

ur\ — —ur2 (at the boundary between Rl and R2) (51) 
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where nm and rifa are the refractive indexes of Rl and R2 regions, re- 
spectively. This reproduces the property (I24|l of our refractive index at 
the boundary. 

In this case, there is no reflection at all for a general angle of incidence 
and polarization between the boundary of Rl and R2. In the same way, 
we can show that there is no reflection between the boundary of R2 and 
R3 (R3 and R4). 

As for the branch cut, when vertical entering (4> = 0) in the case of 
impedance matching [Z\ — Z%), there is no reflection at the branch cut, 
see case (2). Otherwise, there is some reflection, because of the discontin- 
uous behavior of refractive index n and n', see case (3). However, in order 
to achieve perfect invisibility, we can solve this problem by modifying our 
potential U(x) as follows. Due to the striking property of negative refract- 
ing index, the only requirement for potential U(x) in 721 region of the sec- 
ond Riemann sheet is that the light should be bent to the left to the right 
(i.e., the light entering from the first Riemann sheet to Rl region of the 
second Rieman sheet across the branch cut (—2a < x < 2a, y = 0) should 
go through the boundary between Rl and R2 sheet (x = 2a, y > 0)). Our 
potential is the simple example satisfying this condition. However, it is 
easy to see that there are many other suitable potentials as well. There- 
fore, we could slightly modify our potential around the branch cut in Rl 
region to smoothly connecting the potential of the first Rieman sheet to 
the potential of Rl region, keeping the condition that the light should be 
bent from the left to the right. Then, the refractive index can be smoothly 
connected. As a result, there is vanishing reflection around the branch cut 
at all. 

We would like to remark two points. The first remark is that this is 
possible because our proposed design is flexible enough to potentially gen- 
erate a large variety of bounded orbits for light rays using the properties 
of negative refracting indices. The second remark is that this modification 
does not changes the absence of time delay in our device. In other words, 
even if we modify the potential, the time delay is still zero. 

To summarize, we can conclude that perfect invisibility devices in 
isotropic media seem possible. 

7 The size of the invisible space 

In this section, we investigate the maximum size of invisible space. An 
object within the size of this space would be invisible and hidden from 
sight. 

7.1 The radius of invisible space 

Although several local invisible areas may exist, depending on the ray 
light trajectories, we focus on identifying the symmetric size of invisible 
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space characterized by radius R. 

From (|18p . the condition that the r < R in physical space z is invisible 
is as follows. For all point (u, v) of the ray light trajectory in w-space, the 
following equation 

„,2 „,2 

< 1 (52) 



(a 2 /R + R) 2 {a 2 /R-R)'< 

holds. Therefore, the radius of the invisible space is given by the minimum 
of 



r = ^={\/m 2 + v 2 + 4a 2 + \/(u 2 + v' 2 +4a 2 ) 2 — lQa 2 u 2 

-\Jv? +v 2 - Ad 1 + ^[{u 2 + v 2 + 4a 2 ) 2 - 16a 2 u 2 } (53) 

for all points (u, v) of the ray light trajectory 

The possible points for minimizing (153[l are (u,v) = (2a,y ma x) or 
(0,6a), where y max is given by and (JSTJ. 

We will compute these points by considering the following three cases (i), 
(ii) and (iii), respectively. 

(i) In case of (u, v) = (2a, y max ) and (0 < A < J|). 

By substituting u — 2a and v — y ma x into (|53|) . we transform (|53l) into 



ri(j 4) = fe{l + Ji+ . 16 ° a 



Umax y 

Then, substituting (|37() into (J54j) , we obtain 



2 + 2Wl + ^}. (54) 



*w =n a - 4a2{l + v 1 + ^ - y 2 + 2 v 1 + (55) 

In particular, when A — > 

n(A)~^a». (56) 

Thus, when A is small, the radius of invisible space goes to 0. 

(ii)In case of (u,v) = (2a,y max ) and (f^ < A < £■). 
Substituting (|36[) into (|54[) . we obtain 



#m /-, 16a 2 A 2 / „ I 16a 2 A 2 
r*(A) = U {1 + y/l + - y 2 + 2^ 1 + _^-}. (57) 

In particular, this function takes the maximum 

r 2 (^) = (1 + V2 - y2 + 2^2)a ~ 0.21a, (58) 

when A = 

4a 



13 



As an example, we write the value of r 2 for different value of A. When 

A = &T> 

= (2 + \/5 - 2^2 + V5)a ~ 0.12a (59) 
which is smaller than Eq. (|60p . 

(iii) In case of (u, v) = (6a, 0) 

By substituting (u, v) = (6a, 0) into (|53[) . we obtain 

r 3 = (3 - 2^2) = 0.17a. (60) 

By summarizing (i), (ii) and (iii), we obtain the radius of the invisible 
space R as follows: 

f n(A) (for 0<A<|L)) 
i?(A) = I r 2 (A) (for £<A<^)) (61) 
{ r, (for f| < A < £)), 

where n(A), r2(A) and rz are given by (|55|1 . (|57|l and (|60[) . Here, we use 
the fact that {r2(jg-) = ^3)). We show the radius of the invisible space 
(|6ip in Fig. [5] This result indicates that with increasing A, the radius 
of invisible space increases up to 0.17a. As the size of the invisible space 
is proportional to the branch cut a, the proposed invisibility device may 
potentially allow us to hide extensive objects. 



8 Summary 

In this work, we have proposed a new design for a cloaking device made of 
a metamaterial with negative refractive index. Our results indicate that 
by using negative refractive materials, the constraints for designing and 
manufacturing invisibility devices can be more relaxed. In previous works, 
refractive-index profiles on the interior Riemann sheet should guide the 
rays around the branch points and should define close loops or trajectories. 
Therefore, according to classical dynamics [2], the number of potential or 
index profiles was strongly reduced to a few ones: harmonic-oscillator, 
Kepler profile and Maxwell's fish eye. In contrast, our proposed design is 
flexible enough to potentially generate a large variety of bounded orbits 
for light rays using striking properties of negative refractive indices. 

More importantly, our device based on optically isotropic material with 
negatively refracting index represents an example where the time delay is 
zero. Furthermore, due to impedance matching of negatively refracting 
materials, the waves are completely transmitted when crossing different 
material index borders, and reflecting waves are close to zero. These 
results strongly indicate that perfect invisibility with isotropic materials 
is possible. 

We believe that our proposed design may stimulate theoretical and 
practical studies for designing and manufacturing cloaking devices using 
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the properties of negative refraction index. Emerging research based on 
metamaterials has opened up new avenues by exploiting the artificial di- 
electric media. From telecommunications, radar and optical invisibility 
in defense areas to medical imaging and microelectronics in industrial 
sectors, the new emergent electromagnetic features obtained using meta- 
materials can bring exciting scientific progress. Further progress toward 
three-dimensional variations of this scheme is encouraged together with 
experimental work based on cloaking devices for achieving invisibility in 
the visible range of the spectrum. 
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Trajectory of light 
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Figure 1: Trajectory of light rays on the second Riemann sheet in space w. The 
branch is represented by a bold line between -2 to 2 in the real axis. The four areas 
describing the refractive index profile of the cloaking device are shown in the figure. 
Parameters are set to a = 1, E — 4. 
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The maximum radius of invisible space 
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Figure 2: Maximum radius of invisible space. With increasing the potential parameter 
A, the radius of invisible space becomes large up to the value 0.17a. Parameters are 
set to a — 1, E — 4. 
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